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ABSTRACT 

Establishing the reliability with which stellar population parameters can be mea- 
sured is vital to extragalactic astronomy. Galactic GCs provide an excellent medium in 
which to test the consistency of Single Stellar Population (SSP) models as they should 
be our best analogue to a homogeneous (single) stellar population. Here we present 
age, metallicity and a-element abundance measurements for 48 Galactic globular clus- 
ters (GCs) as determined from integrated spectra using Lick indices and SSP models 
from Thomas, Maraston & Korn, Lee & Worthey and Vazdekis et al. . By comparing 
our new measurements to independent determinations we are able to assess the ability 
of these SSPs to derive consistent results - a key requirement before application to 
heterogeneous stellar populations like galaxies. 

We hnd that metallicity determinations are extremely robust, showing good agree- 
ment for all models examined here, including a range of enhancement methods. Ages 
and a-element abundances are accurate for a subset of our models, with the caveat 
that the range of these parameters in Galactic GCs is limited. We are able to show that 
the application of published Lick index response functions to models with fixed abun- 
dance ratios allows us to measure reasonable a-element abundances from a variety of 
models. We also examine the age-metallicity and [a/Fe]-metallicity relations predicted 
by SSP models, and characterise the possible effects of varied model horizontal branch 
morphology on our overall results. 

Key words: globular clusters: general - stars: abundances 



1 INTRODUCTION 

Increased telescope size and improved instrumentation have 
allowed the observation of ever more distant objects. How- 
ever even with these improvements the vast majority of ex- 
tragalactic sources will remain unresolved. The accurate and 
reliable analysis of integrated stellar populations is therefore 
key to our understanding of formation and evolutionary pro- 
cesses in galaxies. Through comparisons of these integrated 
populations with models of homogeneous stellar systems, or 
Single Stellar Populations (SSPs), recent studies have met 
with some success in determining ages and metallicities for 
both galaxies (e.g. Trager et al. 2000; Terlevich & Forbes 
2002; Maraston et al. 2003; Proctor et al. 2004a; Maras- 
ton 2005; Thomas et al. 2005) and extragalactic globular 
clusters (e.g. Forbes et al. 2001; Puzia et al. 2003; Beasley 
et al. 2005; Pierce et al. 2006). 
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Analyses of extragalactic targets have long been depen- 
dent on the accurate modelling of stellar populations. Em- 
pirical approaches to the modelling of integrated light (e.g. 
Spinrad & Taylor 1971) have since given way to more rigor- 
ous models dependent on a knowledge of underlying phys- 
ical processes (i.e. stellar formation and evolution). These 
early analyses primarily made use of broadband colours in 
deriving their age and metallicity measurements. However, 
the limitation of broadband colours is that they are degen- 
erately sensitive to age and metallicity (i.e. old, metal-poor 
and young, metal-rich populations are photometrically iden- 
tical), heavily restricting the accuracy of ages and metallic- 
ities determined using colours alone. 

The addition of spectral indices, in particular Lick index 
absorption features (Burstein et al. 1984; Trager et al. 1998), 
to stellar population models has afforded the much needed 
leverage to break this degeneracy. SSP models including Lick 
indices were first assembled by Worthey (1994), who mod- 
elled 21 Lick indices and sought to identify those features 
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that were particularly age (e.g. Balmer lines) or metallicity 
(e.g. Fe4668, Fe5015, Fe5709 etc.) sensitive and therefore the 
most useful for overcoming the observed degeneracy. More 
recent works (e.g. Maraston 1998; Vazdekis 1999; Bruzual & 
Chariot 2003; Thomas, Maraston & Bender 2003; Thomas, 
Maraston & Korn 2004; Le Borgne et al. 2004; Maraston 
2005; Lee & Worthey 2005) have focused on including more 
indices (i.e. the higher-order Balmer lines H<5 and H7) and 
increasingly complex evolutionary processes (e.g. mass-loss 
and horizontal-branch morphology). 

Along these lines, efforts have also been taken to ac- 
count for known variations in a-element abundance (N, O, 
Mg, Ca, Na, Ne, S, Si, Ti) with respect to Fe-peak elements 
(Cr, Mn, Fe, Co, Ni, Cu, Zn) and their particular effect on 
Lick index measurements. Tripicco & Bell (1995) computed 
the effects of variation in C, N, O, Mg, Fe, Ca, Na, Si, Cr 
and Ti on the 21 Lick indices modelled by Worthey (1994). 
These relative index sensitivities were then used by Trager 
et al. (2000) to modify the SSPs of Worthey (1994), facili- 
tating the measurement of ages, metallicities and a-element 
abundances for a sample of ~40 elliptical galaxies through a 
comparison of H/3, Mgb and <Fe>. Subsequent calculations 
of abundance effects have mimicked the work of Tripicco & 
Bell (1995), adding sensitivity calculations for higher-order 
Balmer lines (e.g. Houdashelt et al. 2002; Korn, Maraston 
& Thomas 2005) and expanding index sensitivities to en- 
compass a broad range of population metallicities (Korn, 
Maraston & Thomas 2005) . 

Studies making use of SSPs generally compare mod- 
elled Lick line-strengths to those measured from integrated 
spectra in order to determine parameters such as age, metal- 
licity and a-element abundance. To have confidence in the 
application of these SSPs to observations, it is important to 
confirm that they can reproduce independently determined 
results from colour-magnitude diagrams (CMDs) and high- 
resolution stellar spectroscopy. Globular clusters (GCs) pro- 
vide a testbed for SSP models as they represent a coeval 
and chemically homogeneous stellar population that should, 
therefore, be analogous to a synthetic single stellar popu- 
lation. The availability of resolved observations for Galactic 
GCs means that accurate ages and metallicities have already 
been determined using CMDs, while a-element abundances 
have been calculated using high-resolution stellar spectra. 

Such an analysis was carried out by Proctor 
et al. (2004b), who fit a sample of 24 Galactic GC spec- 
tra from Cohen et al. (1998) and Puzia et al. (2002) to SSP 
models from Vazdekis (1999), Bruzual & Chariot (2003) and 
Thomas, Maraston & Bender (2003) using a multi-index \ 2 - 
minimization technique, as opposed to 2-dimensional fits e.g. 
Trager et al. (2000). They found that it was possible to re- 
cover the known age, metallicity and a-element abundance 
to within ~0.1 dex. However their small sample size (20 in- 
dividual GCs) and the relative lack of independent age and 
[a/Fe] determinations mean a reliable statistical comparison 
to literature trends was not possible. 

In this study we expand the work of Proctor 
et al. (2004b), comparing high signal-to-noise (S/N~100) 
spectra of 42 Galactic GCs from Puzia et al. (2002) 
and Schiavon et al. (2005) to recent SSP models from 
Thomas, Maraston & Korn (2004), Lee & Worthey (2005) 
and Vazdekis et al. (2007). To these models we apply a- 
element enhancement calculations from both Houdashelt 



et al. (2002) and Korn, Maraston & Thomas (2005). 
We are then able to determine ages, metallicities and a- 
enhancements using the integrated spectra. A comparison 
between our determined values and those from CMD stud- 
ies (e.g. De Angeli et al. 2005) and resolved stellar spectra 
(e.g. Pritzl et al. 2005) gives an indication as to the relia- 
bility of parameters derived solely from integrated spectral 
analysis. 

Section [3 contains a brief summary of each of the mod- 
els, describing their specifics and the means by which they 
have been calibrated. In Section|3]we detail the Galactic GC 
spectral data used in this work. This section also contains a 
discussion of the a-element enhancement models used and 
an outline of their application. Section|4]describes the multi- 
index x 2 minimization technique that we use for fitting spec- 
tra to the SSP models. In Section [S] we discuss comparisons 
of parameters we derive from SSPs and literature values. 
Conclusions are given in Section [S] 



2 STELLAR POPULATION MODELS 

Models from Thomas, Maraston & Korn (2004; TMK04), 
Lee & Worthey (2005; LW05) and Vazdekis et al. (2007; 
V07) have been chosen for application to the Galactic GC 
data. TMK04 and LW05 models are both computed us- 
ing the Worthey et al. (1994) fitting functions and provide 
Lick/IDS system index values for a range of ages and metal- 
licities. Vazdekis et al. provide their models as SEDs, from 
which we then measure Lick indices. We have decided not 
to include the commonly used models of Vazdekis (1999) or 
Bruzual & Chariot (2003) as they have been studied in a 
similar fashion by Proctor et al. (2004b). 

Below is a summary of the essential parameters for 
each of the SSP models selected for study in this paper. 

Thomas, Maraston &: Korn (2004; TMK04): These 
are based on previous work by Thomas, Maraston & 
Bender (2003; TMB03). Models cover the metallicity range 
— 2.25<[Z/H]<0.65 with ages from 1 to 15Gyrs and are 
based on isochrones from Cassisi, Castellani & Castellani 
(1997), Bono et al. (1997) and Salasnich et al. (2000). 
TMK04 include horizontal branch effects, providing em- 
pirically calibrated Balmer lines modelled for both red 
and blue horizontal branch morphologies using the stellar 
mass loss parameter r\ (Reimers 1975). Variations in 
abundance ratios are tabulated using updated response 
functions that include the higher-order Balmer lines and a 
metallicity dependence as calculated by Korn, Maraston & 
Thomas (2005). These models cover all 25 Lick indices in a 
wavelength range of AA4000-6500 A. Data can be found at 
http : // www . dsg . port . ac . uk/ "thomasd/ . 

Lee & Worthey (2005; LW05): The Lee & Worthey mod- 
els cover a metallicity range of — 2.5<[Fe/H]<0.3 and an age 
range of 1 to 12 Gyrs. Recent Y 2 isochrones (Yi et al. 2001; 
Kim et al. 2002) are adopted along with post-red giant 
evloutionary tracks from Yi et al. (1997). An additional scal- 
ing factor X) is used to account for stellar mass loss and aids in 
matching observed horizontal-branch morphology in Galac- 
tic GCs. SSPs include alpha enhancements of [a/Fe]=0.0,0.3 
and 0.6 applied at super-solar metallicities using updated 
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response functions from Houdashelt et al. (2002). At 
sub-solar metallicities, a-element abundances are super- 
solar and reflect the local abundance-ratio pattern, which 
includes some metallicity dependence. The SSPs model 
25 Lick indices from H<5a to TiC>2, and can be found at 
http : //astro . wsu . edu/hclee/wpRGB_all_Lick_2005. 

Vazdekis et al. (2007; V07): Models from Vazdekis 
et al. are based on the previous models of Vazdekis (1999) 
and Vazdekis et al. (2003). These models are presented as 
SEDs and cover a metallicity range of — 2.3<[Z/H]<0.2 
and ages from 0.1 to 17.5 Gyrs using Padova group 
isochrones from Girardi et al. (2000). These models are 
derived using the recent MILES spectral library (Sanchez- 
Blazquez et al. 2006). Non-solar abundance ratios are 
not accounted for in these models, so SSPs represent the 
local abundance pattern. SED models are available at 
http : / / www . iac . es/galeria/ vazdekis/. 



2.1 SSP Model Calibrations 

For the SSP models outlined above, calibrations have been 
carried out in order to verify the accuracy of their index pre- 
dictions. This is a key step in the construction of these SSP 
models, as the results obtained from their use on extragalac- 
tic sources are generally blind (i.e. there are no corroborat- 
ing methods like CMDs or resolved spectroscopy available) . 
For their calibrations, both TMK04 and LW05 make use of 
the P02 Galactic GC data as measured using the Worthey 
et al. (1994) index definitions. The P02 observations were 
taken with specific care given to their luminosity sampling 
in order to obtain accurate cluster spectra with account of 
stochastic effects. This careful sampling means that spectra 
are representative of the total cluster population and there- 
fore ideal for the calibration of Lick index models. 

Calibrations of the TMK04 models have been well doc- 
umented in Maraston et al. (2003), TMB03 and TMK04, 
which involve assuming an old GC age (12 Gyrs) and com- 
paring measured GC indices to SSP predictions. In Maras- 
ton et al. (2003), these comparisons are carried out using 
index-index comparison with <Fe> (iron-sensitive indices) 
or Mg b (all other indices) and metallicity comparison with 
CMD [Fe/H] determinations (their Figs. 1,7-11). Additional 
evaluations of the higher-order Balmer lines are carried out 
through US, H7 vs. [MgFe] index comparisons (their Fig. 
13). In all of these Maraston et al. (2003) note that the GC 
data of P02 lie as expected in relation to their SSPs, pre- 
dicting metallicities and a-element abundances consistent 
with those from CMD and resolved spectral studies. These 
results are reiterated in TMB03, and the Balmer lines are re- 
calibrated in TMK04, with the same good agreement being 
found. 

Lee & Worthey (2005) perform similar calibrations for 
their SSP models, comparing Lick index measurements of 
the CBR98 and P02 datasets to their SSP models. This is 
done through a comparison of Lick indices to [Fe/H] as pre- 
dicted by their SSPs, using metallicities from Harris (1996) 
for the GC data. They, like KMT04, find good agreement 
between their SSPs and GCs and note that their models 
require no zero-point offset to match the GC data. 



At the time of writing, information regarding the cali- 
bration of V07 models was unavailable. 



2.2 Non-solar abundance ratios 

An important consideration in fitting our sample of Galactic 
GCs is the handling of non-solar abundance ratios. It is well 
known that GCs exhibit elemental abundances that differ 
from those measured in the Sun (Pilachowski et al. 1983; 
Gratton 1987), and tabulated response functions have al- 
lowed for these variations to be included in the SSP models. 
Tripicco & Bell (1995; TB95) modelled response functions 
for the Lick/IDS index system, providing fractional index 
variations for 21 Lick indices with respect to 10 elements 
(C, N, O, Mg, Ca, Na, Si, Ti, Cr and Fe) in three different 
stellar types (cool dwarf, main-sequence turnoff dwarf and 
cool giant). The TB95 fractional responses were calculated 
by doubling each element, Xi, in turn ([Xi/Fe]=+0.3) and 
measuring the resultant effects on each index. While TB95 
calculations were carried out using a 5 Gyr old isochrone, 
adjusting the relative contribution of their three modelled 
stellar types allows the construction of stellar populations 
with a range of ages, metallicities and a-element ratios. 

Work presented by Houdashelt et al. (2002; H02) has 
sought to update the response functions of TB95. They have 
used recent, updated line lists to improve upon the original 
TB95 calculations and include the higher-order Balmer lines 
(R8 a ,f and H7a,f) and TiO not modelled by TB95. H02 
have also varied the method by which carbon enhancement 
is calculated. Rather than double carbon (+0.3 dex) they 
have chosen to enhance carbon by only +0.15 dex, seeking to 
avoid modelling discrepancies that arise as C/O approaches 
1. 

Additional work by Korn, Maraston & Thomas (2005; 
KMT05) has tested some of the simplifications made by 
TB95. KMT05 confirm the validity of performing all calcu- 
lations using a 5 Gyr isochrone by comparing to the results 
of calculations made using a 1 Gyr isochrone. In this com- 
parison they find only small deviations between the two re- 
sults, ~1 percent for G4300 and Fe4348 and significantly less 
for all other indices. In addition, they include calculations 
for high-order Balmer lines (H<5a,f, H7a,f), TiO and have 
added metallicity dependence to their fractional responses, 
recalculating the same tables as TB95 for six different metal- 
licities from -2.25<[Z/H]<+0.67. 

In this work we will be using both the KMT05 and 
H02 fractional sensitivities applied to a variety of models 
using the methods described in Trager et al. (2000; TOO) and 
TMB03. For specific details of this enhancement application 
and calibration, we refer the reader to Appendix A. 

For clarity, SSPs that are used as originally published 
will be referred to by their respective references (i.e. TMK04, 
LW05 and V07), while models that we have altered through 
the use of the H02 and KMT05 index response functions will 
be referred to by a combination of the model and enhance- 
ment calculation reference (e.g. TMK+H02, LW+KMT05, 
V+H02 etc.). 
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Figure 1. Mean deviations of each index for each model set. Indices clipped in the iterative process described in Section [4] are not 
included. Circles and open squares represent fits to H02 and KMT05 enhanced models respectively. For V07 based models H/3 was found 
to be particularly deviant and excluded from all fits to these models. It is shown here for comparison only. 



Table 1. Number of times each index is clipped for each model 
set and enhancement method. Results for H02 enhanced models 
are shown, with results for KMT05 enhanced models shown in 
parentheses. 



Lick 
Index 


LW05 


TMK04 


V07 


a SA (A) 





(0) 





(0) 





(0) 


n 5F (A) 





(0) 





(0) 


3 


(3) 


G4300 (A) 


(i 


(5) 


9 


(8) 


8 


(18) 


H^a (A) 


1 


(0) 





(0) 





(1) 


H-yF (A) 





(0) 


1 


(1) 





(4) 


Fe4383 (A) 


2 


(1) 


1 


(1) 





(1) 


Ca4455 (A) 


48 


(44) 


44 


(36) 





(1) 


Fo4531 (A) 





(0) 





(0) 





(0) 


C4668 (A) 


33 


(38) 


28 


(25) 


34 


(39) 


H (A) 


1 


(2) 


4 


(1) 






Mgi (mag) 


7 


(8) 


5 


(6) 


4 


(2) 


Mg 2 (mag) 


1 


(0) 


2 


(12) 


1 


(3) 


Mg 6 (A) 


8 


0) 


14 


(23) 


1 


(ID 


Fo5270 (A) 


2 


(1) 


2 


(3) 


2 


(0) 


Fe5335 (A) 





(1) 


3 


(3) 


21 


(14) 


Fe5406 (A) 


1 


(0) 





(0) 


2 


(2) 


Fo5709 (A) 





(1) 


1 


(1) 


1 


(2) 


Fo5782 (A) 


20 


(24) 


27 


(20) 


3 


(0) 


TiOi (mag) 


1 


(1) 





(0) 


1 


(1) 


TiQ 2 (mag) 


6 


(4) 


1 


(1) 


2 


(2) 



Index measurements for the S05 spectra were carried 
out using index definitions from Worthey et al. (1994), 
Worthey & Ottaviani (1997) and Trager et al. (1998). Prior 
to index measurement, spectra were broadened to the Lick 
system resolution (~8-llA) using a wavelength dependent 
Gaussian broadening kernel based on the IDS resolution de- 
scription given in Worthey & Ottaviani (1997). Measure- 
ments for 25 indices were produced, however unreliable 
fluxes around 4546 A and 5050 A from CCD defects or 
sky subtraction errors resulted in deviant measurements for 
the Fe4531 and Fe5015 indices (see S05 for details). 

The second source of Galactic GC data is from the study 
of Puzia et al. (2002; P02). P02 provide long-slit spectra for 
12 Galactic GCs in the wavelength range AA3400-7300 A. 
Observations were carried out using the European Southern 
Observatory (ESO) 1.52-m telescope on La Silla with the 
Boiler & Chivens Spectrograph with a spectral resolution of 
~6.7A per pixel. Lick line-strengths are given for 25 indices 
measured and calibrated using both the Trager et al. (1998) 
and Worthey et al. (1994) index definitions. These data will 
be discussed further in Section [3]2] See P02 for more details 
regarding observations and line-strength measurements. 



3 GALACTIC GC SPECTRAL DATA 

The Galactic GC spectra used in this study are taken from 
two different sources. The first data are from Schiavon 
et al. (2005; S05) who provide spectra for 41 GCs chosen to 
represent a range of parameters (e.g. age, metallicity, Galac- 
tocentric distance, etc.). These spectra were obtained at the 
Cerro Tololo Inter- American Observatory (CTIO) Blanco 
4-m telescope with the Ritchey- Chretien spectrograph and 
cover a wavelength range of AA3350-6430 A at a resolution 
of ~3.1 A per pixel for central wavelengths. This allows mea- 
surement of all 25 Lick indices from H5 a to Ti02- For addi- 
tional details regarding observations, see S05. 



3.1 Multiple observations 

The combined dataset used for this analysis includes spec- 
tra for 42 unique GCs with 75 observations in total. P02 
contains no duplicate spectra, however S05 include multiple 
observations and aperture extractions for several GCs. In or- 
der to assemble a more coherent sample, analyses have been 
limited to a single observation for each GC in each study. 
In the case of multiple observations, we have selected to use 
those that are best fit (i.e. most indices fit with the lowest 
X 2 ) for the majority of SSP models. For GCs with multiple 
aperture extractions (NGC 6284, NGC 6342, NGC 6441, 
NGC 6528, NGC 6624 and NGC 7078) we have used the 
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extraction that includes a wider spatial region than just 
the FWHM of the slit profile as these, generally, give the 
best fit. We find deviations in parameters derived from fits 
across multiple observations to be small (±0.015, ±0.031 
and ±0.038 in log age, [E/Fe] and [Z/H] respectively), and 
consequently exclude them from our analysis has little effect 
on our final results. 

3.2 Calibration to the Lick/IDS system 

The general method for calibrating observations to the 
Lick/IDS system involves obtaining spectra of stars in the 
Lick standard library and using these to calculate line- 
strength offsets. Such calibrations have been carried out for 
the P02 dataset, however S05 observed Jones library (Jones 
1999) standard stars which, owing to the slightly limited 
spectral coverage of the Jones library (AA3820-5410 A with 
a gap from 4500 A to 4780 A), only allow for the calibra- 
tion of at most 17 Lick indices. In an effort to include as 
many indices as possible in our SSP model fitting, we have 
instead chosen to calibrate the S05 dataset using the 11 GCs 
it shares with P02. 

Puzia et al. (2002) provide Lick indices measured using 
both the Worthey et al. (1994; W94) and Trager et al. (1998; 
T98) passband definitions. The differences between W94 and 
T98 index definitions are the result of refinements to the 
wavelength solution of the original Lick/IDS library spec- 
tra and constitute 1.25 A to 1.75 A shifts in index defini- 
tions. Central indices (H/3 to Fe5406) were unaffected by 
this adjustment as their original definitions were calculated 
using more finely calibrated template spectra (Worthey 
et al. 1994). 

Because of the index adjustment made by T98, their in- 
dex definitions are the most appropriate for use on properly 
wavelength calibrated data; index and pseudo-continuum 
passbands will fall on the correct spectral features. For this 
reason the "correct" index definitions to use for the P02 
data are from T98, however concerns have been raised with 
regards to the calibration of these index measurements to 
the Lick/IDS system. Specifically, P02's W94 measurements 
were calibrated using index values published by the Lick 
group, while their T98 indices were calibrated using indices 
re-measured from the published spectra (D. Thomas & C. 
Maraston 2006, priv. comm.). Further examination of the 
P02 data shows large offsets between their W94 and T98 
measurements, even for indices whose passband definitions 
remain the same between W94 and T98. 

In light of these inconsistencies we have performed cal- 
ibrations of the S05 data using both the W94 and T98 in- 
dex definitions. Indices were measured on the S05 spectra 
using both the W94 and T98 passband definitions and cali- 
brated using the corresponding data from P02. This method 
of Lick calibration introduces a greater uncertainty in our 
calculations than if Lick standard stars were used. We have 
therefore adjusted our index errors accordingly, including 
the rms about the mean offsets for the common GCs and 
the rms quoted in P02 for their own calibration to the Lick 
system (from their Tables 3 and Dl) in our overall error es- 
timates. The final errors we adopt for this calibration are 
shown in Appendix B, Table iBll 

For the remainder of this work we to show fits to the 
W94 calibrated data for TMK04 and LW05 based models 



to avoid the uncertainties in P02's calibration of their T98 
data, discussed above. All relevant figures have been repro- 
duced using the T98 calibrated data for comparison and are 
shown in Appendix [Cl 

3.3 Vazdekis 2007 Models 

Data fit to the V07 models do not require the same Lick/IDS 
calibrations as data fit to the TMK or LW based models, 
since V07 use a well calibrated stellar library. With this in 
mind, the S05 data were broadened using the same wave- 
length dependent Gaussian discussed above, and indices 
were measured using the T98 index definitions. No addi- 
tional calibration was performed. 

In order to fit the P02 data to the V07 models, we 
have used the coefficients given for the Lick calibration of 
their data measured using T98 definitions (P02's Table 3) 
to de-calibrate the indices given in their Table Cl. This re- 
sults in indices measured using T98 passband definitions on 
smoothed, flux calibrated spectra. 

In Appendix [B] we show index-index comparisons for 
the common GCs in P02 and S05 for the three difference 
calibrations discussed above. 



4 GALACTIC GC FITS USING SSP MODELS 

Having detailed the models and data used, we now turn to a 
discussion of our fitting technique. We have chosen to adopt 
the x 2 -fitting method discussed by Proctor et al. (2004b), 
involving the simultaneous ^-minimization of as many in- 
dices as possible in order to maximize use of the available 
data and break the age-metallicity degeneracy. This tech- 
nique has been used previously to determine ages and metal- 
licities of GCs (Proctor et al. 2004b; Beasley et al. 2005; 
Pierce et al. 2005; Pierce et al. 2006), and has been shown 
to produce more robust results than most individual index 
comparisons. The method involves the ^-minimization G f 
measured spectral indices to a grid of SSP indices corre- 
sponding to different metallicities, ages and a-element abun- 
dances. Indices that show significant deviations (~3a) from 
the best fit may be removed and the fits recalculated. This 
process can be continued until no more deviant indices are 
present and a stable fit is established. The fits produced by 
this method are robust against single deviant indices and 
calibration errors, and allow for the reliable identification of 
trends across multiple data sets. 

The iterative fitting and clipping of indices involved 
in this multi-index technique makes easy the identification 
and omission of indices that are deviant for a majority of 
the GC spectra. The NaD index, for example, is known 
to suffer heavily from interstellar absorption, and so ex- 
hibits large variations across the data sets when fit; we have 
therefore excluded this index from all fits. As is commonly 
the case in GCs, we find that residuals to best fits of the 
CN and Ca4227 indices follow a pattern suggestive of ni- 
trogen enrichment (CNi and CN2 show positive residuals, 
while Ca4227 shows a negative residual; e.g. TMB03, Proc- 
tor et al. 2004b). Rather than fit nitrogen as an independent 
parameter (e.g. TMB03) , we have simply excluded these in- 
dices from the fitting procedure. 

For most of the S05 spectra, measurements of the 
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Figure 2. SSP grids of H/3 vs. [MgFe] for TMK04, LW05 and V07. Circles (blue) and squares (red) represent S05 and P02 data 
respectively. Grid lines cover metallicities from [Z/H]=0.0 (solar) down to —2.20 in steps of 0.25 dex. Age lines are 1, 1.5, 2, 3, 5, 8 and 
12Gyrs for LW05, with 15 Gyr lines included for TMK04 and 17Gyr lines for V07. Arrows indicate the direction of increasing age or 
metallicity. The average error is shown in the upper right corner (see Section l3.2|l . 



Fe4531 and Fe5015 indices were inhibited by "deviant 
fluxes" in their index bands, attributed to poorly subtracted 
sky lines or CCD defects. Both Fe4531 and Fe5015 have been 
excluded from all fits to the S05 spectra (See S05 for more 
details). 

As discussed by Proctor et al. (2004b) the Fe5015 
index showed large deviations between the P02 and Co- 
hen, Blakeslee & Ryzhov (1998) data sets, which are pos- 
sible symptoms of the inconsistencies in conversion to the 
Lick/IDS system mentioned in Section[3] Again, since there 
is no absolute way to account for these deviations, the 
Fe5015 index was excluded from fits to the P02 data. Taking 
all of these effects into account, we are left to conduct fits 
using 20 indices for the P02 data, and 19 indices for S05. 

Table [1] gives the details of the fits, showing the to- 
tal number of times that each index was clipped for a given 
model set. Ca4455 was found to be very deviant in the LW05 
and TMK04 models and can be interpreted as a problem in 
calibrations to the Lick/IDS system, as fits to the V07 mod- 
els (un-calibrated data) do not show the same deviations. 
For all 6 model sets the C4668 index is particularly deviant, 
perhaps due in part to its extreme carbon sensitivity (the 
adopted enhancement pattern leaves carbon solar scaled). 
However the improved fit of this index to V07 models that 
some part of this offset could also be due to calibration (as 
with Ca4455). We found the H/3 index to be particularly 
aberrant in fits to V07 models, and so this index has been 
removed from all fits to V07 (both P02 and S05 data). 

In Figure [T] we show the mean deviations for our best 
fits of each index, for each set of models. Fits to models using 
H02 enhancements (TMK+H02, LW+H02 and V+H02) are 
shown as filled circles and fits to KMT05 enhanced models 
(TMK04, LW+KMT05 and V+KMT05) are show as open 
squares. Enhancement methods give qualitatively the same 
fits (to within a fraction of the errors) for a given model set. 
Comparing the quality of fit between LW05 or TMK04 and 
V07, one sees the benefit of fitting to models based on a 
well calibrated stellar library, in this case MILES, reflected 
in the reduced deviations seen for most indices. 

In all subsequent figures, fits to TMK+KMT05 mod- 
els are shown as open symbols plotted behind the parame- 
ters derived using TMK04 models, which are shown as filled 



symbols. This is done for comparison only, and analyses are 
carried using TMK04 models results. 



5 COMPARISON OF SSP DERIVED 
PARAMETERS WITH LITERATURE 

We have measured ages, metallicities and a-element 
abundances for Galactic GCs using several stellar popu- 
lation models. In doing this, we are able to compare the 
SSP derived parameters to those determined using other 
methods (i.e. CMDs or resolved stellar spectroscopy) and 
assess the validity of SSP determinations. Establishing the 
reliability of these SSP predictions is vital as these models 
are frequently applied to extragalactic sources (galaxies and 
GCs) for which alternative age, metallicity and a-element 
abundance determinations are not available. It is important 
to note that we have therefore conducted fits to age, metal- 
licity and a-element abundance simultaneously, rather than 
assuming an old age (e.g. Maraston et al. 2003; TMB03), 
in order to duplicate the way in which these models are 
frequently used for extragalactic sources. 

Fig. [2] shows model grids of H/3 plotted against [MgFe] 
for TMK04, LW05 and V07. GC data from S05 and P02 are 
overplotted. From Fig. [2] both TMK04 and LW05 models fit 
the data reasonably well, with the largest deviations gener- 
ally at intermediate metallicities (-1.0<[Z/H]<-0.5) where 
horizontal branch morphology becomes increasingly uncer- 
tain (see Section r5.5|l . The comparison of V07 models to data 
is less encouraging, as the data lie well below the grids at 
nearly all metallicities. This is consistent with our findings 
from the x 2 fits, namely that the H/3 index is particularly 
deviant when compared to other indices. 

At the low metallicity end, we see that models differ in 
their coverage of the observed data. However as it is at these 
metallicities that stellar libraries become extremely sparse, 
this variability is not unexpected. In all cases data are con- 
sistent (within errors) with the models, which is important 
to our fitting procedure as stable, accurate fits are difficult 
to obtain for GCs whose index values fall outside the range 
of the models. 

Having briefly looked at the base models, we now turn 
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Figure 3. [Fe/H] from Harris et al. (1996) plotted against SSP derived metallicities. Numbers in the upper left corner represent the 
mean offset from the one-to-one line (dashed) and the cr rms scatter about that offset, error bars signify a 1<t deviation on our SSP fits 
and ± 0.1 dex for the Harris [Fe/H] values. 



to a discussion of the ages, metallicities and a-element abun- 
dances derived using these SSP models. 



5.1 Metallicity 

A direct comparison between TMK04, LW05 and V07 mod- 
els is not straight-forward, as each have handled metallicity 
in a slightly different way. TMK04 present their models in 
terms of [Z/H], the total metallicity, and have accounted 
for the local stellar pattern in their models (e.g. Wheeler 
et al. 1989), and so the values of [Z/H] and [E/Fe]jmeasured 
using their SSPs can be used without any adjustment. In 
contrast, LW05 models are supplied as a function of [Fe/H]. 
As LW05 models do not, as published, account for varying 
a-element ratios at [Fe/H]<0, these models carry with them 
an implicit enhancement, [a/Fe]i oca i, equivalent to the local 
stellar abundance pattern (i.e. [a/Fe]=0.3 for [Fe/H]<-1.0; 
[a/Fe] decreasing from 0.3 to 0.0 as [Fe/H] increases from 
-1.0 to solar; [a/Fe]=0.0 for [Fe/H]>0.0). This pattern must 
be accounted for in our measurements, in addition to the 
enhancement [E/Fe]ssp that we measure from our own en- 
hancement calculations (see Appendix A). For LW05 mod- 
els, [Z/H] is then calculated with Equation Al using our 
measured [Fe/H] and [a/Fe]=[a-/Fe]i 0ca i + [E/Fe] ssp • 

The V07 models provide yet another variation, being 



1 [E/Fe] is used to describe the measured enhancement, as it 
represents an enhancement of N, O, Mg, Ca, Na, Si and Ti as 
opposed to any single element. 



published as a function of [Z/H], but do not include vary- 
ing a-element ratio calculations. Since the [Z/H] we measure 
is already includes the afore mentioned local stellar abun- 
dance pattern, and our additional enhancement calculations 
are applied at constant total metallicity, [Z/H], we can then 
calculate [Fe/H] in the same way as LW05, using Equation 
Al and [a/Fe] = [a/Fe] locat + [E/Fe] ssp. 

Compounding these modelling differences is the uncer- 
tainty of the Zinn & West (1984; hereafter ZW84) GC metal- 
licity scale used in Harris (1996). While ZW84 is gener- 
ally quoted as [Fe/H], it is based on measurements made 
by Cohen et al. (1983) using the average of the Mg triplet 
(~5175A), 5270 A and 5206 A Fe blends. The output of a 
particular SSP, then, is somewhat coloured by assumptions 
made as to what Harris (1996) actually represents (be it iron 
abundance, [Fe/H], or overall metallicity [Z/H]). This ambi- 
guity is magnified by evolutionary tracks and stellar libraries 
which may or may not have made additional assumptions as 
to the nature of ZW84 [Fe/H] values. 

In light of these ambiguities, in Fig. [3] we show [Fe/H] 
measurements from Harris (1996) plotted against the most 
closely related metallicity indicator from each of the mod- 
els. For TMK04 and V07 models, this is [Z/H] ssp, however 
for LW05 this is [Fe/H] ssp- This represents a fundamen- 
tal difference in what the models are measuring, and should 
be kept in mind when these models are applied to spectra. 
This difference is likely due to the several factors mentioned 
above, however it does not prevent a qualitative compar- 
ison of these models. In fact all model variants shown in 
Fig. [3] find metallicities (either [Z/H] s ST for TMK04 and 
V07 or [Fe/H] ssp for LW05) that are in good agreement 
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Figure 4. Zinn & West (1984) scale ages from De Angeli et al. (2005) plotted against SSP derived ages, the dotted line in this case 
representing the oldest age in each model set. Symbols represent P02 (squares) and S05 (circles). Numbers in the upper left corner 
represent the mean offset from the one-to-one line (dashed) and the cr rms scatter about that offset, error bars signify a la deviation on 
our SSP fits. 



with the CMD metallicities from Harris (1996), regardless 
of enhancement method. Looking more closely at the offsets 
and scatters for each panel in Fig. [3] (A and a rms respec- 
tively in the upper left corners), the TMK04 and TMK+H02 
models give the tightest relations (i.e. lowest oy ms ), albeit 
with a slightly larger offset from the one-to-one line than 
LW05 based models. V07 models do not seem to follow the 
one-to-one line as closely as the other 4 models, generally 
showing larger offsets and scatters than either of the other 
two models. 



5.2 Age 

While the range of Galactic GCs ages is quite small 
(~2Gyrs), a comparison of SSP age predictions to liter- 
ature is still useful in evaluating their reliability. In Fig. 
[4] we show such a comparison, with CMD determined 
ages from De Angeli et al. (2005; hereafter D05) plotted 
against SSP ages. Average ages are 10.74±1.84, 9.38±1.82 
and 11.70±3.60 Gyrs for TMK+H02, LW+H02 and V+H02 
models respectively; for models using KMT05 enhancement, 
mean ages are 10.78±1.63, 9.60±1.79 and 11.47±3.16 Gyrs. 

Models based on V07 SSPs do a comparatively poor job 
of GC age prediction, finding mean offsets from D05 ages and 
scatters about these offsets (A and a rms in the upper left 
corner of each panel in Fig. [4} significantly larger than ei- 
ther TMK04 or LW05 based models. Modelling uncertainties 
in the age-sensitive Balmer lines initially seemed the likely 
culprit for these large deviations, however an examination of 
the H<5 and H7 indices did not show a significant offset (i.e. 



as is observed in the H/3 index for these models; discussed in 
Section [3J. As an additional test to this, fits to V+H02 and 
V+KMT05 were conducted with all Balmer lines omitted 
(H5a,f, H7a,f and H/3), however no significant change was 
observed, i.e. ages were still found to be abnormally high 
with large scatter. As these age deviations are present in 
many indices (i.e. more than just the age sensitive Balmer 
lines), it seems that either data need some additional cali- 
bration to be properly fit to V07 models, or an additional 
calibration of the models themselves is needed. 

TMK04 and LW05 based models do a good job of repro- 
ducing D05 ages, both finding reasonably small mean offsets 
(A<lGyr). LW+KMT05 models do the best quantitative 
job of reproducing the CMD ages of D05, giving both the 
smallest mean offset (A=-0.057) and scatter (oy ms = 1.544), 
however there are several important caveats to this age anal- 
ysis. 

Firstly, differing upper age limits for each of the models 
(15, 12 and 17 Gyrs for TMK04, LW05 and V07 respectively; 
dotted lines in Fig. 2]) likely play some role in the apparent 
agreement or disagreement of SSP ages with literature val- 
ues. Most notably, fits to LW05 based models find several 
GCs with ages equivalent to the upper limit of the models, 
whereas in fits to TMK04 derived models, all GC ages are fit 
as opposed to being assigned the maximum available value. 

It should also be noted that differences in modelling, es- 
pecially evolutionary tracks, can affect the age comparisons 
shown in Fig. 3] While Lick indices should not be affected 
by the particular set of evolutionary tracks used (Maraston 
et al. 2003), varied handling of a-element abundance ratios 
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will influence the agreement between SSP and CMD derived 
ages. Both De Angeli et al. (2005) and LW05 use a-enhanced 
isochrones (Cassisi et al. 2004 and Kim et al. (2002) re- 
specively), and so the good agreement in their predicted 
ages could be a result of this. Conversly, both TMK04 and 
V07 models use solar scaled evolutionary tracks, which have 
been shown to produce slightly older age estimates than 
a-enhanced isochrones of a similar metallicity (Salasnich 
et al. 2000). 

Additionaly the CMD ages of De Angeli et al. (2005) are 
subject to uncertainties in their absolute calibration, being 
similarly based upon model isochrones. The ages shown in 
Fig. 2] should therefore be viewed as measuring the relative 
agreement of two difference methods of age measurement, 
CMD vs. spectral, rather than a comparison of absolute 
ages. 



5.3 a-Element abundance 

The ability of SSP models to accurately measure a-element 
abundances is of great interest as it can give an indication 
of formation timescales in galaxies. Measurements of en- 
hancement using SSPs, [E/Fe]ssp, are shown in Fig. ^plot- 
ted against [Mg/Fe], [(Ca+Ti)/Fe0 and [(Mg+Ca+Ti)/Fe] 
from Pritzl et al. (2005). The agreement between Pritzl 
et al.'s [Mg/Fe] and [E/Fe]ssp is poor, likely owing to the 
inclusion of additional elements in the SSP enhancement 
[E/Fe]ssp- Perhaps not surprisingly, both [(Ca+Ti)/Fe] and 
[(Mg+Ca+Ti)/Fe] from Pritzl et al. (2005) relate more 
closely to [E/Fe] SS p- From Fig. [5] we see that KMT04 and 
LW+H02 are able to best reproduce the enhancement val- 
ues of Pritzl et al. (2005), with [(Mg+Ca+Ti)/Fe] being the 
best fit while TMK+H02, V+KMT05 and V+H02 appear 
to underpredict GC element abundances. 

We find mean [E/Fe] values of 0.12±0.08, 0.28±0.11 
and 0.24±0.06 for TMK+H02, LW+H02 and V+H02 re- 
spectively. Models using KMT05 enhancement calcula- 
tions produce higher [E/Fe] SS p=0.28±0.13, 0.37±0.12 and 
0.22±012. All models, with the exception of TMK+H02, 
produce mean [E/Fe] ssp values consistent with literature 
findings of a constant a-element abundance of [E/Fe]~0.3 
(e.g. Gratton et al. 2004). 

To further examine model enhancement predictions, in 
Fig. HJ] [E/Fe]ssp is plotted against the same SSP metallic- 
ity indicators as Fig. [3] ([Z/H] for TMK04 and V07, [Fe/H] 
for LW05). The shaded region in Fig. [^represents the range 
covered by field star data from Pritzl et al. (2004; their Fig. 
4) . This abundance trend with metallicitjjf] is generally at- 
tributed to the increased influence of Type la SNe at later 
times when higher metallicity stars formed, and so it is not 
surprising that GCs formed at a similar epoch (i.e. metallic- 
ity) are found to follow this same a-enhancement pattern. At 
higher metallicities, i.e. [Fe/H]>-0.5, evidence for GCs ex- 
hibiting the same a-enhancement "down-turn" as field stars 
is less certain, however this is largely due to the small num- 
ber of high-metallicity GCs relative to lower metallicities. 



2 [a/Fe] from Pritzl et al. (2005) 

3 This trend is generally shown as a comparison with [Fe/H] , 
however here we have plotted [E/Fejggp against the indepen- 
dently determined metallicity measure for each model. 
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Figure 5. SSP derived values of [E/Fe] plotted against high res- 
olution element abundances from Pritzl et al. (2005). Models and 
symbols are the same as in Figure [3] Error bars signify a la devi- 
ation on our SSP fits and ±0.1 dex for the high-resolution [a/Fc] 
values. 



With regards to the SSPs fit here, TMK04, LW+H02 
and V+H02 most closely match the observed field star 
abundance pattern. LW+KMT05 models seem to over- 
predict enhancements at all metallicites (as evident from 
the mean [E/Fe] SSP =0.37), while TMK+H02 greatly under- 
predict enhancement at low metallicities and are incon- 
sistent with the field star pattern. All KMT05 enhanced 
models show some deviation at low metallicities, either 
over-predicting (i.e. TMK04 and LW+KMT05) or under- 
predicting (i.e. V+KMT05) enhancement. The variation 
seen in Fig. [6] between TMK04 and LW+KMT05, which 
find higher [E/Fe]ssp at low metallicities, and V+KMT05, 
which find lower [E/Fejssp at low metallicities is a result of 
the V07 models being offset from the data (as seen in Fig. 
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Figure 6. Comparison of fitting results for TMK04, LW05 and V07 based models. Symbols are the same as in previous figures. The 
shaded region shows the range of abundances observed in field stars from Pritzl et al. (2005). GCs fit at the minimum SSP metallicities 
are not shown. 



[2|. These deviations, however, are accompanied by increas- 
ing error distributions and so are roughly consistent with 
the bulk of the data. This will be discussed in more detail 
below. 



5.3.1 Houdashe.lt et al. (2002) vs. Korn, Maraston & 
Thomas (2005) Enhancement 

As previously mentioned, one marked difference between the 
KMT05 and H02 enhanced model sets is the tendency for 
KMT05 enhanced models to show odd enhancement behav- 
ior at low metallicities, be it the increased [E/Fc] values at 
low metallicities in the TMK04 or LW+KMT05 models or 
the abnormally low [E/Fe] values in V+KMT05 models. 

The primary difference between the KMT05 and H02 
enhancement calculations is the inclusion of metallicity de- 
pendent index sensitivities by KMT05. In their calculations, 
KMT05 found that indices at low metallicities are relatively 
insensitive to variations in a-element abundance. This insen- 
sitivity means that all of the SSP grid lines "pinch" together 
at low metallicities, resulting in indices that may be only 
slightly enhanced in line-strength relative to [E/Fe]=0.0 be- 
ing measured as having elevated enhancement (as is the case 
in both TMK04 and LW+KMT05 models) and increased er- 
ror distributions (e.g. V+KMT05). 

To test that the primary difference between H02 and 
KMT05 enhanced is, in fact, the added metallicity depen- 
dence rather than an overall shift in enhancement calcu- 
lation between H02 and KMT05 we have constructed two 
new sets of models (using LW05 and TMK04) using the 
[Z/H]=0.0 index sensitivities from KMT05 applied at all 
metallicities, making them comparable to H02 enhanced 
models. When we compare ages, metallicities and [E/Fe] val- 
ues between these models and their H02 enhanced counter- 
parts (TMK+H02 and LW+H02) we find that differences 
are of order an interpolation step (±0.025 dex in log(age) 
and metallicity, ±0.03 dex in [E/Fe]). Deviations between 
KMT05 and H02 models, then, are almost entirely due to 
the metallicity dependence added by KMT05. 



5.4 Age-Metallicity Relation 

The age-metallicity relation (AMR) for Galactic GCs is 
well established and shows that Galactic GCs are generally 
old (e.g. Salaris & Weiss 2002; Beasley et al. 2005; Puzia 
et al. 2005; De Angeli et al. 2005). De Angeli et al. (2005) 
have most recently examined the AMR of Galactic GCs us- 
ing HST imaging and found that very low metallicity GCs 
([Fc/H]<-1.4) are old (~llGyrs) with very low scatter in 
their ages (0.06 Gyrs), however at intermediate metallicities 
GC ages show considerably more variety in their ages, rang- 
ing from 7.5 to 11 Gyrs. 

While the precision of SSP age measurements is not 
fine enough to delineate between slight variations in age 
(our general uncertainty is ~2 Gyrs), it is of interest to test 
whether SSP models can reproduce the overall trend of uni- 
form, old GC ages. In Fig. [7] we plot the AMR as derived 
from each of the SSP models. 

Models based on TMK04 produce AMRs that reveal 
an odd trend of increasing age towards higher metallicity. 
It should be noted, however, that age determinations at 
the lowest metallicities are highly uncertain due to possible 
variations in the horizontal branch morphology. As a result, 
while there is a suggestion of a positive AMR slope, the data 
are consistent with a uniformly old GC system. TMK04 es- 
pecially does an excellent job of reproducing a generally old 
GC population, finding only one GC younger than 8 Gyrs. 

LW±H02 and LW±KMT05 model fits show GCs con- 
sistently old at metallicities [Fe/H]<-1, however they have 
a high metallicity "tail" towards younger ages. CMD age 
determinations for these highest metallicity GCs, NGC 6528 
and NGC 6553, suggest that their actual ages are old (11- 
13 Gyrs, Zoccali et al. 2001; Feltzing & Johnson 2002) and 
consistent with the rest of the Galactic GC system. Compar- 
ing LW05 based models with TMK04, we note that TMK04 
does not appear to have the same problem at high metal- 
licities, finding ages for NGC 6528 and NGC 6553 consistent 
with CMD determinations. Errant age measures from LW05 
based models could be indicative of issues with this par- 
ticular model set at higher metallicities, but we are unable 
to comment in more detail on the quality of fits at higher 
metallicites due to a lack of data points. At metallicities 
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Figure 7. [Fe/H] vs. age as derived from our SSP fitting. Symbols and models are the same as in previous figures. The dotted line 
represents the maximum age for a particular SSP. 



[Fe/H]<-0.4 the LW+H02 model is in excellent agreement 
with CMD based AMRs (e.g. De Angeli et al. 2005). 

V07 based models have difficulty in producing an 
AMR consistent with the known GC AMR (e.g. De An- 
geli et al. 2005). This is almost certainly due to the age 
determination problems discussed earlier (see Section [5.2[) . 
which makes any useful AMR determination nearly impos- 
sible. Again, the large age variations seen when fitting to 
these models (^lOGyrs) may be alleviated via calibration 
to the MILES stellar library. 

5.5 Horizontal Branch Morphology 

Balmer line indices (RSa,f, H7a,f and H/3) are particularly 
sensitive to the presence of hot stars, becoming weaker as 
temperatures decrease. This lends to their use as age indi- 
cators as the decrease in main-sequence turnoff luminosity, 
and therefore temperature, associated with an aging stel- 
lar population is echoed strongly in the measured Balmer 
line strengths. However as older populations are considered 
(>10Gyr), the increased presence of hot horizontal branch 
(HB) stars causes Balmer line strengths to increase, lead- 
ing to an age degeneracy at low metallicities, with very old 
stellar populations appearing young. 

The modelling of these HB morphologies is particularly 
difficult, as the interplay of contributing effects (e.g. mass- 
loss, metallicity, dynamical effects etc.) is not known well 
enough to be modelled in detail (i.e. based purely on theory). 
Modelling varying HB morphologies, then, has been done 
primarily via prescriptive methods. Maraston & Thomas 
(2000) find that a mass-loss prescription is able to reproduce 
the strong Balmer lines found in old elliptical galaxy popu- 



lations, as well as the trends of increasing H/3 line strengths 
in low-metallicity Galactic GCs. 

Of the SSP models used here, only TMK04 allow for a 
variation of HB mophology in their models, supplying two 
sets of empirically calibrated Balmer line indices (H<5a,f, 
Hya.f and H/3), one each for Red and Blue HB morphol- 
ogy (see Maraston & Thomas 2000 for details). In Fig. \E\ 
we show a comparison of the BHB (solid lines) and RHB 
(dashed lines) grids supplied by TMK04. This figure illus- 
trates the need for a consideration of variable HB morphol- 
ogy, especially at intermediate metallicities where hot hori- 
zontal branch stars begin to cause large variations in Balmer 
line strengths. 

In an effort to better quantify the HB effects in our re- 
sults, we have performed a second set of fits to the RHB 
models of TMK04 using the same techniques described in 
Section. U In Fig. [5] we show the difference between age, 
metallicity and a-element abundance derived using BHB 
and RHB models plotted against the GC horizontal branch 
ratio (HBR) from Harris 1996 and Zoccali et al. (2000). Since 
the HBR is based purely on numbers of stars in a given 
branclQ it is an excellent, independent means of determin- 
ing HB morphology. 

The first thing to note in Fig. is that GCs with 
HBRs>0.9 are almost entirely blue, and so have been ex- 
cluded from the statistics shown in the bottom left of each 
panel. As the RHB modelling in TMK04 is limited to Balmer 
line indices, it is not surprising that the changes seen in 
[Z/H] and [E/Fe] as a result of this modelling are small. In 
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Figure 8. TMK04 SSP grids for blue (solid lines) and red 
(dashed lines) horizontal branch morphologies. P02 (circles) and 
S05 (squares) GC data overplotted. 



particular, the offsets for both are dominated by the scatter. 
Ages are most strongly affected by varying HB morphology 
in the TMK04 SSP models, however even mean offset for 
these (~0.11) is of order the error for our age determina- 
tions (~0.1). 

While the results of this comparison do show that HB 
morphology is important to individual indices (e.g. H/3 in 
Fig. [8} , it does no appear that it greatly affects the proper- 
ties that we derive using the multi-index fitting technique. 
Most importantly, the offsets that we see as a result of vary- 
ing the HB modelling are not large enough to significantly 
change the results of our analyses. 



6 CONCLUSIONS 

Using integrated spectra of Galactic GCs, we have 
tested SSP model predictions of age, metallicity and a- 
element abundances using Lick indices. The multi-index \ 2 - 
minimization technique used has allowed us to measure GC 
stellar population parameters consistent with published val- 
ues, even in situations where data are poorly fit in a single 
index-index space. 

Although Galactic GCs cover a wide range of metal- 
licities (-2.28<[Fe/H]<-0.04 in our sample), our GC sam- 
ple only represent a limited parameter space in age 
and a-element abundance (the mean age of our GCs 
is 10.23±0.83 Gyrs from De Angeli et al. 2005, and 
[a/Fe]=0.30±0.09dex from Pritzl et al. 2005). To expand 
this work, SSP model parameter space should be explored 
further by obtaining integrated spectra for Galactic GCs 
covering a broader range of ages (e.g. Terzan 7 or Pal 12 with 
intermediate ages) and a-element adundances (e.g. Rup 106 
or Pal 12 with sub-solar a-element abundances). 

We find metallicity to be the most robust parameter, 
showing almost no sensitivity to the different enhancement 
calculations. We note that differences in the construction of 
models have lead to a fundamental difference in the metal- 
licities they predict. In particular, [Z/H] measurements from 
TMK04 and V07 models are mostly closely related to CMD 
[Fe/H] from Harris et al. (1996). However for LW05 models, 
SSP predictions for [Fe/H] most closely match the metallic- 
ities of Harris et al. (1996). This does not affect our conclu- 
sions, however is an important caveat to keep in mind when 
appling these models. 

Age determinations using TMK04 and LW05 based 
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Figure 9. Difference in metallicity, enhancement and age be- 
tween red and blue horizontal branch models plotted against hor- 
izontal branch ratio (HBR) from Harris (1996). Symbols are the 
same as previous figures. Offsets and scatters are shown in the 
lower left corners for all GCs with HBRs<0.9. 



models are very reliable for the old GCs of the Milky Way. 
However, V07 models have difficulty in recovering reliable 
age measurements, giving GCs that are too old and with a 
very large scatter. 

Of the models tested here, only TMK04 models pro- 
vide measurements of a-element abundances at all metal- 
licities. However, we have shown that with a relatively sim- 
plistic application of Lick index sensitivity calculations from 
either Korn, Maraston & Thomas (2005) or Houdashelt 
et al. (2002) we are able to recover reasonable a-element 
abundances (i.e. consistent with the literature) from each 
SSP model. 

We have shown that HB morphology is an important 
consideration, as it can dramatically affect the age sensitive 
Balmer line indicies. However, in spite of these Balmer line 
variations, we find changes in SSP determined ages as a re- 
sult of varying HB morphology (i.e. blue or red HB models) 
are relatively small for the majority of GCs when using the 
multi-index \ 2 fitting method. Determinations of metallicity 
and a-element enhancement are relatively robust to changes 
in HB morphology, a key result for extragalactic GC studies 
where direct determinations of HB morphology are unavail- 
able. 
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APPENDIX A: DETAILS OF ENHANCEMENT 
APPLICATION 

Here we provide some more details of our ad hoc enhance- 
ment applications involving the calculations of H02 and 
TMK05. In summary, these methods involve the adjustment 
of ct-element abundances at a fixed metallicity. Using the re- 
lation 



[Z/H] = [Fe/H] + A[a/Fe] 



(Al) 



which, applied at constant metallicity leads to the following 
equation from (TOO): 



A = 



A[Fe/H] 
' A[a/Fe] 



(A2) 



where A varies depending on which elements are selected 
as enhanced or depressed. In this work we consider N, O, 
Mg, Ca, Na, Si and Ti as a-group elements, while Cr and 
Fe are depressed. We leave C unchanged (i.e. solar-scaled). 
This enhancement scheme mimics the work of TMK04, al- 
lowing us to compare our TMK+KMT05 models to pub- 
lished TMK04 models that use these same index sensitivity 
calculations. Using Equation IA2I we find A=0.934 when C 
is excluded from the enhanced group and solar-scaled. Fol- 
lowing TMB03, the enhancement applied using the following 
equation: 



= Iss P Y[eMRo.s(i)] (AlXt]/ °- 3) 



(A3) 



where the quantity flo.3(ijj refers to the index change 
resultant from increasing the abundance of the i th element 
by 0.3 dex (adopting notation from TOO). This equation 
is arrived at by assuming In I is a linear function of [Xi] 
and using a Taylor expansion to approximate the effects of 
abundance ratio changes (see TMB03 for details). 



Al Negative indices 

One of the problems that arises when applying abundance 
ratio adjustments is that, in some cases, Lick indices have 
negative values. In applying the expansion of In I (Equation 
IA3[) however, it is implicit that the index values cannot be 
negative, and in fact asymptotically approach zero at low 
abundances (TMB03; TOO). KMT05 handle this problem by 
applying their computed index response directly to the flux 
of their absorption lines, yielding a positive result. Lack- 
ing the means to apply corrections in this manner, we have 
adopted the method used by TMB03, in which the lowest 
value of a particular index at a given age is taken as the 
zero point, and all other indices are scaled to reflect this 
zero point shift. In the case of the high-order Balmer lines 
this occurs at the highest metallicity, while for most other 
Lick indices it occurs at the lowest. For C4668 and NaD, 
the low points do not occur at either end of the metallic- 
ity scale, indicating an inflection point in the indices at an 
intermediate metallicity. For these indices, values from the 
local minima are adopted. The adjustment, using notation 
of TMB03, is defined as follows: 



n 1 1 1 1 i 1 
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Figure Al. Our adopted contributions of each stellar type with 
age as approximated using Fig. 41 of Worthey (1994). This allows 
for the luminosity of older populations to be increasingly affected 
by stars occupying the giant phase of evolution, and proved im- 
portant in matching the behaviour of the TMK04 models. 



8 = lion 



\hou 



(A4) 



This process is summarized in equation 9 from TMB03: 

(A[X«]/0.3) 



" ( 1 dl V 



(A5) 



where 5 is now the zero-point index adjustment, and Io is the 
absolute index value taken from the index response tables. 
The key difference between TMB03's and our application 
of this zero-point adjustment is that TMB03 adjust each 
evolutionary phase individually. Here we are left to apply 
this adjustment to the integrated index measurements. 



A2 Additional Considerations 

In calculating our fractional index sensitivities, we have in- 
cluded several additional effects not discussed in Sections 
EUlandEI] Firstly, the fractional responses of KMT05 have 
been adopted strictly in a differential sense, using the SSP 
index values (I ssp ) in place of their Io values. While shift- 
ing from Io to I ssp should generally improve accuracy, for 
indices in which the lowest index value is very small (i.e. 
Issp — S < 0.05) the calculated fractional responses are over- 
estimated, since i?o.3(i) c< (I SS p — <5) _1 . To account for this 
we have applied a correction to Equation IA4I such that 



8 — Ilo-w \Ilou 



(A6) 



1 81 

Io d[Xi] 



0.3 



This adjustment increases all indices relative to the zero 
point, and should have little effect on the applied fractional 
enhancements except in the cases of very low index values. 

Secondly, additional considerations were included to ac- 
count for variations in the flux contribution of a given evolu- 
tionary type with respect to both age and wavelength. This 
was done using estimated values from Fig. 41 from Worthey 
et al. (1994). The results of our approximations are shown 
in Figs. IA1I and IA2I The effect of this added age and wave- 
length dependence is generally 5-10 percent, however can be 
as much as 50 percent for the bluest indices in the youngest 
populations. We stress that these are approximations only, 
made to include some handling of varying flux contributions 
for different populations. 
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Table Al. Columns 2 and 3 give the a Tma between [E/Fe]=0.3 
TMK04 and TMK+KMT05 models (see Section \M} in two dif- 
ferent metallicity ranges. In most cases the errors in our enhance- 
ment application are well within our Lick calibration errors and, 
therefore, have little effect on the derived parameters. 
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Figure A2. Our adopted contributions of each stellar type as 
they varies with wavelength. These show the same trends as noted 
by Worthey (1994), that for younger populations the equivalence 
point for giant and main sequence stars shifts to higher wave- 
lengths. 
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Since the published TMK04 models make use of the 
KMT05 enhancement calculations, our goal in applying en- 
hancements was to mimic the TMK04 models with our own 
TMK+KMT05 calculations. To this end, we have compared 
the TMK+KMT05 models to the TMK04 models for twice 
solar ([E/Fe]=0.3) grids. The results of this comparison are 
shown in Columns 2 & 3 ([Z/H]<0 and [Z/H]<0 respec- 
tively) of Table |A"T1 We find the agreement to be acceptable, 
particularly at the lower metallicities where the majority of 
our GCs lie. In nearly all cases deviations we find are less 
than our adopted errors for conversion to the Lick/IDS sys- 
tem and should, therefore, have negligible effects on parame- 
ters (age, metallicity and a-enhancement) derived from fits. 
The C4668 index is an exception to this, showing significant 
deviations in the higher metallicity regime. This is indicative 
of the high sensitivity of the C4668 index to small variations 
in a-element abundance. 
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Table Bl. Summary of errors associated with the calibration of our GC data. Adopted Lick/IDS calibration errors for S05 measured 
and calibrated using Worthey et al. (1994) and Trager et al. (1998) index definitions are shown in Columns 2 and 3 respectively. 
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APPENDIX B: INDEX-INDEX CALIBRATIONS 

Here we show the results of calibrating the S05 data to the P02 data. Table IB 1 1 shows the adopted error we have associated 
with these calibrations. Figures [Bll IB2I and IB3I show index-index comparisons for the 11 GCs common between P02 and S05 
for each of the three different calibrations. 



APPENDIX C: DATA CALIBRATED USING TRAGER et al. (1998) INDEX DEFINITIONS 

As discussed in Section [3T2J P02 provide data measured and calibrated using both the W94 and T98 index definitions, however 
in the text we use only the W94 calibrated data. Here we reproduce Figs. [3H3 using the P02 and S05 data measured and 
calibrated using T98 index definitions. 
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Figure Bl. A comparison of the common GCs between S05 and P02 using the Worthey et al. (1994) index definitions for measurement 
and calibration. Solid lines are a one-to-one correlation, with dashed lines representing our adopted Lick/IDS conversion error as shown 
in Column 2 of Table |BT1 Indices not included in fits to the SSP models (see Sect. [4} are shown in grey. 
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Figure B2. Identical to Fig, IBll but with indices measured and calibrated using the Trager et al. (1998) index definitions. Solid lines 
are a one-to-one correlation, with dashed lines representing our adopted Lick/IDS conversion error as shown in Column 3 of Table IBTI 
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Figure B3. Identical to Fig, IB1I but with indices measured using the Trager et al. (1998) index definitions, but without any additional 
Lick calibration applied. Solid lines are a one-to-one correlation, with dashed lines representing our adopted Lick/IDS conversion error 
as shown in Column 4 of Table IB1I 
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Figure CI. [Fc/H] from Harris et al. (1996) plotted against SSP derived metallicities. Numbers in the upper left corner represent the 
mean offset from the one-to-one line (dashed) and the cr rms scatter about that offset, error bars signify a 1<t deviation on our SSP fits 
and ± 0.1 dex for the Harris [Fe/H] values. 
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Figure C2. Zinn & West (1984) scale ages from De Angeli et al. (2005) plotted against SSP derived ages, the dotted line in this case 
represents the oldest age in each model set. Symbols represent P02 (squares) and S05 (circles). Numbers in the upper left corner represent 
the mean offset from the one-to-one line (dashed) and the oy ms scatter about that offset, error bars signify a lcr deviation on our SSP 
fits. 
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Figure C3. SSP derived values of [E/Fe] plotted against high resolution element abundances from Pritzl et al. (2005). Models and 
symbols are the same as in Figure ICll Error bars signify a lcr deviation on our SSP fits and ± 0.1 dex for the high-resolution [a/Fe] 
values. 
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Figure C4. Comparison of fitting results for TMK04, LW05 and V07 based models. Symbols are the same as in previous figures. The 
dashed line shows the assumed local abundance pattern in stars. GCs fit at the minimum SSP metallicities are not shown. 
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Figure C5. [Fe/H] vs. age as derived from our SSP fitting. Symbols and models are the same as in previous figures. The dotted line 
represents the maximum age for a particular SSP. 



